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Abstract 

We consider tunneling transitions between states separated by an energy barrier in 
a simple field theoretical model. We analyse the case of soliton creation induced by 
collisions of a few highly energetic particles. We present semiclassical, but otherwise 
first principle, study of this process at all energies of colliding particles. We find that 
direct tunneling to the final state occurs at energies below the critical value E c , which is 
slightly higher than the barrier height. Tunneling probability grows with energy in this 
regime. Above the critical energy, the tunneling mechanism is different. The transition 
proceeds through creation of a state close to the top of the potential barrier (sphaleron) 
and its subsequent decay. At certain limiting energy E\ tunneling probability ceases to 
grow. At higher energies the dominant mechanism of transition becomes the release of 
energy excess E — E\ by the emission of a few particles and then tunneling at effectively 
lower energy E = E\ via the limiting semiclassical configuration. The latter belongs to 
a class of "real-time instantons" , semiclassical solutions saturating the inclusive proba- 
bility of tunneling from initial states with given number of particles. We conclude that 
the process of collision-induced tunneling is exponentially suppressed at all energies. 



1 Introduction 

Tunneling processes are inherent in many field theoretical models, well-known examples 
being false vacuum decay in scalar theories pQ and topology-changing transitions in gauge 
theories |2] . To some extent tunneling in field theory is similar to that in quantum mechanics. 
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A relevant energy scale for tunneling is set by the minimum height of the energy barrier 
separating initial and final states. It is equal to the energy of a critical bubble Ij and 
sphaleron [3] in the two examples above. Classical transition is energetically forbidden 
at energies below the barrier height. In the weak coupling regime semiclassical technique 
is applicable for the description of tunneling, and the rate of transition is exponentially 
suppressed. On the other hand, one might expect the suppression to vanish once large enough 
amount of energy, lifting the system above the potential barrier, is injected. This expectation 
turns out to be correct in some situations, such as transitions at finite temperature (see 
Ref . j3] and references therein) , at non-zero fermion density [S] , or in the presence of heavy 
particles in the initial state jHllZl- There is, however, a situation which seems to provide a 
notable exception from this rule. This is the case of tunneling induced by collision of two 
highly energetic particles. 

A prospect of observation of non-perturbative transitions in particle collisions [HI Ej pro- 
voked extensive discussion in literature (see [TUl HU Ej for reviews) . Semiclassical results 
|12[ show, however, that this process of induced tunneling remains exponentially 

suppressed even when the energy of colliding particles exceeds considerably the height of 
the potential barrier. Furthermore, investigations of toy models ITo] and unitarity argu- 
ments JdESEI] indicate that the process of induced tunneling should remain exponentially 
suppressed even as the energy of collision tends to infinity. Nevertheless, no direct calcula- 
tion of the suppression exponent at very large energies of colliding particles has been carried 
out so far in any model, not to speak about the limit of infinite energy. 

In this paper we consider a simple field theoretical model allowing for the semiclassical 
calculation of the suppression exponent of induced tunneling at all energies of colliding 
particles. The model describes free scalar field </>(£, x) living in (1+1) dimensions on a half- 
line x > 0, with interactions localized at the boundary point x = 0. The action of the model 
is 

oo 

S = ^JdtJdx [{d^f - m 2 2 ] - 4 J dt [1 - cos(#(t, 0))] . (1) 



The second term represents boundary interaction with the characteristic energy scale \i. The 
bulk mass m is introduced as an infrared regulator, and it is assumed to be small compared 
to the scale \x. In the main body of the paper we take the limit m — > 0, as the small mass 
turns out to be irrelevant for the semiclassical evaluation of the suppression exponent. The 
only role it plays is to fix the free particle basis for the asymptotic states. 

The massless version of the model with different infrared regularizations has been 
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considered in several contexts. In condensed matter physics it has been used to describe 
transport in quantum wires (201 an d Josephson chains with defects |2T], the role of the 
regulator in those cases being played by the length of wire or chain. In [22] , a special value 
of the coupling constant g = \[2~tx was considered and the model was shown to be exactly 
solvable in that case. Finally, in [2S] the model (0) with m = 0, obtained as a limit of a model 
with sine-Gordon interaction in the bulk, was argued to be integrable at any value of g. 

The model (JTJ) possesses a number of static solutions determined by the value of the field 
at x = 0. Up to corrections of order m/fi, these solutions have the following form: the 
value of the field at the boundary is "pinned down" to one of the minima of the boundary 
potential, 0^ = 2irn/g, n = 1, 2, . . . , while in the bulk the field slowly falls off due to the 
presence of the small mass, 

These solutions are localized near the boundary, x = 0, so it is natural to call them "boundary 
solitons" . The masses of the solitons, 

M« = ^-m , (2) 
9 

are relatively small, as they are proportional to the bulk mass m. 

In this paper we study the process of creation of the first soliton (<f)^ = 2n/g) in a collision 
of a particle (or a bunch of particles) with the boundary 3 ; the energy of the particle(s) is 
assumed to be much larger than the soliton mass. To start with, let us consider the classical 
analog of this process, which is the creation of the soliton in the collision of a classical wave 
packet with the boundary. It is easy to see that the classical process is possible only if the 
energy of the wave packet exceeds some threshold energy Eg, which is much higher than 
the soliton mass (|2"|. Indeed, the boundary value <f>(t, 0) of the relevant classical solution 

changes from to 2ir/g during the process. The classical solution passes the maximum of 

( 

the boundary potential <fi — ir/g at some moment of time, and thus the total energy of 
this solution is larger than the boundary energy at the maximum, 

Es = %. (3) 

g 2 

3 If the bulk mass m is equal to zero (and the model is somehow regularized in infrared) the spatially 
homogeneous configurations 4>^(x) = (jx^ are classical vacua of the model. Our results then apply to 
vacuum-to-vacuum transitions induced by collisions of highly energetic particle(s) with the boundary. All 
formulae in the main text remain valid. Yet another situation where our analysis applies literally, is induced 
false vacuum decay in a model with m = 0, and degeneracy between the classical vacua <^ n )(x) = 
slightly lifted. 
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One concludes that any state containing boundary soliton is separated from the vacuum by 
a potential barrier. It is straightforward to find static configuration "sitting" on top of the 
barrier, which we call "sphaleron" following Ref. [3J. This is an unstable static solution of 
the classical field equations representing the saddle point of the static energy functional. In 
our model it has the same exponential form as the soliton but with the boundary value on 
top of the boundary potential, 

9 

The energy of the sphaleron is given by Eq. (jSJ), again up to corrections of order m/yU. 

We see that the soliton production in collisions of particle(s) with the boundary is clas- 
sically forbidden and hence exponentially suppressed at energies smaller than the sphaleron 
energy. The question is what happens when energy grows. In this paper we study this 
question by applying semiclassical methods. The semiclassical approximation is justified by 
the following observation. After the rescaling of the field, (f) —>■ <p/g, the coupling constant g 
enters the action only through the overall multiplicative factor 1/g 2 . Therefore, g 2 plays the 
role of the Planck constant h, and the weak coupling limit corresponds to the semiclassical 
situation. This is the case we consider in this paper. 

The techniques we use are explained in more detail in subsequent sections. Here we briefly 
summarize our results. We calculate semiclassically the probability V of the soliton produc- 
tion in collision of a particle with the boundary. In the leading semiclassical approximation 
it has the exponential form, 

and we concentrate on the calculation of the suppression exponent F . Our results for the de- 
pendence of the suppression exponent on the energy of incoming particle are shown in Fig. ^ 
There are three energy regions corresponding to three physically different mechanisms of the 
process. If the collision energy is lower than the critical value E c « l-2.Es (region I of the 
figure), tunneling occurs in conventional way with the relevant semiclassical configurations 
ending up directly in the soliton sector. This regime can be called "direct tunneling". We 
study this region of energies in Sec. 01 and find the analytic result for the suppression of such 
transitions: 

E < E c . (4) 

Formula (@J), if continued to the energies higher than the critical energy E c , would show that 
the transitions become unsuppressed at energy i^Eg. However, this formula is incorrect at 



F(E) = 4vrln 



E 



4 



k 1 

1 \ 

1 


1 1 1 

II 


III 


1 

1 1 1 


» — - 



0.5 1 E c 1.5 2 2.5 Ey 3 

E/E s 



Figure 1: Dependence of the suppression exponent F on the collision energy E, the latter 
measured in units of the sphaleron energy Eg = 2fi/g 2 . The tunneling process is driven by 
three physically different mechanisms in the parts I, II, III of the graph. 

E > E c , as the solutions describing direct tunneling cease to exist at energies higher than 
the critical one. 

Semiclassical configurations with energies above E c are studied in Sec. 0] using a combi- 
nation of analytic and numerical methods. We find that their properties are different from 
the ones below the critical energy, and the tunneling mechanism they describe is entirely 
different. Instead of tunneling directly to the other side of the energy barrier, the system 
jumps on its top, thus creating the sphaleron which then decays producing the soliton in the 
final state. The second stage of the process, decay of the sphaleron into the soliton, proceeds 
with the probability of order one. Still, the transitions remain exponentially suppressed 
due to considerable rearrangement the system has to undergo during the first stage of the 
process, i.e. formation of the sphaleron. The tunneling mechanism outlined here has been 
observed recently in quantum mechanics [2E] and gauge theory Transitions of this type 
correspond to the region II of Fig. [TJ they extend up to the energy E\ « 2.7 Eg- 

The solution describing the transition with the energy E\ belongs to a novel class of 
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semiclassical solutions which we call "real-time instantons" . The latter saturate the inclusive 
probability of tunneling from initial states with given number of particles. We consider the 
real-time instantons in Sec. 14.31 

The energy E\ is the best energy for the transition. Further increase of energy above 
E[ (region III of Fig. ^) does not lead to any gain in the transition probability, and the 
suppression exponent stays constant all the way up to E = oo. We show in Sec. 14.41 that the 
system emits the energy excess (E — E{) in the form of a few highly energetic particles, and 
then undergoes the transition of the second type (on top of the energy barrier) with the energy 
Ei. Emission of one or several particles does not change the suppression exponent, affecting 
only the pre-exponential factor. In this way we obtain that F(E > E{) = 10.27. Our results 
confirm the conjecture [21] that the suppression exponent of collision-induced tunneling is 
frozen at high energies. The tunneling mechanism at E > Ei, i.e. emission of energy excess in 
the form of a few particles and tunneling via the most favourable semiclassical configuration, 
coincides with that conjectured in Ref. |15|- 

Section |5] contains concluding remarks. 



2 T/6>-problem 
2.1 General formalism 

A difficulty one encounters in the semiclassical description of collision-induced tunneling 
is that the initial state of the process is not semiclassical. To this end, the method of 
calculation should be appropriately adjusted. In this work we adopt the method proposed 
in |2Sj- Namely, let us consider the inclusive probability of tunneling from multiparticle 
states with energy E and number of particles N: 

V(E,N) = J2\(f\SPEpN^\ 2 , (5) 

where <S is the S'-matrix while Pe and P/v are projectors onto states with given energy and 
number of particles. The states \i) and |/) are perturbative excitations above the vacuum 
and the soliton, respectively. The function (j5J) can be calculated with the use of semiclassical 
methods, provided that the energy and initial number of particles are semiclassically large, 
E = E/g 2 , N = N/g 2 . The result has the exponential form, 



V(E,N) oce 



-F(E,N)/g 2 
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It is clear that the multiparticle probability V(E, N) provides an upper bound on the tunnel- 
ing probability induced by collisions of fewer particles. Indeed, any initial few-particle state 
can be promoted to a multiparticle one by adding an appropriate number of "spectator" 
particles which do not affect the tunneling process. Thus, the exponent F(E) of the soliton 
production in the collision of one particle with the boundary, is larger than the inclusive 
suppression exponent F(E,N). Furthermore, one conjectures [25 that the small-particle 
limit of the multiparticle exponent coincides with the exponent of tunneling induced by one 
or a few particles, 

F(E) = lim F(E,N) . (6) 

The relation (JBJ) has been checked in several orders of perturbation theory around the instan- 
ton in gauge theory j2E], and also non-perturbatively in quantum mechanics of two degrees 
of freedom [23 I2E|- We use Eq. (JHJ) throughout this paper to obtain F(E) from the mul- 
tiparticle exponent F(E,N). To simplify notations, we omit tilde over the rescaled energy 
and number of particles below. 

Transitions with fixed initial number of particles N are described by solutions of the so- 
called T/9 boundary value problem j2H|. The latter is formulated on the contour ABCD in 
complex time plane shown in Fig. |21 Namely, the configurations describing tunneling should 
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Figure 2: Contour in complex time on which the boundary value problem is formulated, 
satisfy the classical field equations in the internal points of the contour, 

(d^ -dl + m 2 )0 = 0, x > , (7a) 
<9 x (/> = sin 0, x = 0. (7b) 

The Euclidean part BC of the contour may be interpreted as representing the evolution of 
the field under the barrier, "duration" T of this evolution is a parameter of solution. Field 
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equations (|7a|) . ()7b|) are supplemented by initial and final boundary conditions in the parts 
A and D of the contour, respectively. Namely, the field <ft(t,x) should be real as t — > +00; it 
describes the final state of the semiclassical evolution, 

Im — > as t — > +00 . (7c) 

The field in the initial state is linear about the vacuum = 0, 

1 f dk 



4>{t' + iT,x) 



V— ► — OO 



and the boundary conditions in the part A of the contour relate positive and negative fre- 
quency components of the solution, 

fk = e~ 9 9k ■ (7e) 

The boundary condition (|7e|) can be understood as follows. In the limit 9 — > +00 it coincides 
with the Feynman boundary condition and thus corresponds to the initial state with semi- 
classically small number of particles. Finite 9 corresponds to the initial state with non-zero 
N, from which tunneling occurs in the least suppressed way. 

Given the values of T and 9, one finds a complex solution <fi(t, x; T, 9) satisfying equa- 
tions (|7aj) — (ffejl . The energy and initial number of particles for this solution are given by 
the familiar formulae, 



E = J dk uj k fk9*k > ( 8a ) 
N = [dk f k g* k . (8b) 



Alternatively, they can be determined by differentiating the action functional evaluated on 
the solution, with respect to the parameters T and 9, 

E = — lmS(T,9), N = 2—lmS(T,9), (9) 



S= \ I dt I dx [-0<9^0 - (<9 X 0) 2 - m 2 2 ] -/I I dt (1 - cos0(t, 0)) (10) 







where 

S = 

2 _ 

is the action of the model 4 , integrated by parts and calculated along the contour ABCD. 
The suppression exponent of the process is given by the Legendre transform of the action 
functional, 

F(E, N) = 2lmS — N9 — 2ET . (11) 
4 Hereafter we use the rescaled action which does not contain the coupling constant g. 
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Below we refer to the problem (J7J) as "T/0-problem", and use the term "#-instanton" for 
the relevant semiclassical solution. 

Several remarks are in order. First, the boundary value problem (jZJ) does not guarantee 
that its solutions interpolate between states with and without the soliton. In the next 
subsection we describe additional requirements that should be imposed to ensure that the 
solutions are relevant for tunneling. 

Second, one observes that the solution <p(t, x) can be analytically continued to the com- 
plex time plane, and in this way the contour ABCD may be deformed without affecting the 
integral (fTTj) for the suppression exponent. The only thing one should worry about while 
deforming the contour is to avoid the singularities of the solution, which are shown schemati- 
cally by the double lines in Fig. |21 Below it will be convenient not to be attached to a contour 
of any particular form. Instead, we look for the solution <fi(t, x) satisfying Eqs. (JTaJ), (ffbj) 
in the entire complex time plane, with boundary conditions (|Tc|) and (fTej) imposed in the 
asymptotic regions D and A of the complex plane. When using this approach, one should 
make sure, however, that the asymptotic regions A and D can be connected by a contour 
avoiding the singularities of solution. 

Finally, let us discuss two particular limits of the boundary value problem ((7|). The 
first one corresponds to 9 = 0. Solutions obtained in this limit are periodic in Euclidean 
time, they are called "periodic instantons" (2S1- These solutions determine the inclusive 
probability of tunneling from states with given energy and arbitrary number of particles, 

V p (E) = J2\(f\SPE^)\ 2 • 

*,/ 

Periodic instantons have been extensively studied in literature [3U], in particular, they can 
be used to determine the rate of transitions at finite temperature. 

Below we will encounter solutions of the T/#-problem obtained in the limit T = 0. To 
the best of our knowledge, such solutions have never been considered before. We call them 
"real-time instantons", as the contour ABCD does not contain a Euclidean part in this case. 
These solutions depend on a single auxiliary parameter 9, i.e. a single free physical parameter 
of the incoming state, which we choose to be the number of particles N. The energy is thus 
a function of N for this set of solutions. The value of the suppression exponent F rt (N) 
evaluated on the real-time instantons is non-zero exclusively due to the fact that these field 
configurations are complex- valued. At first sight, the latter property may seem to contradict 
the condition (J7c|) requiring reality of the solution in the region D of the real axis. However, 
this is not the case: the condition (J7c|) is imposed only in the asymptotic future, and it does 
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not prevent the solution from being complex-valued at finite times. Let us discuss this point 
in more detail. Assume that at large (but finite) time the solution is linearized about some 
static configuration, 

0(t,x) = O (x) + 50i(t,x) + i502(t,x) , 

where both 50i and 502 are real. The condition (J7c|) implies that 502 — ¥ at t — > +oo. This 
amounts to requiring that 502 describes the evolution along the negative mode of the static 
solution 0o( x ), so that 502 oc e _const '*. Thus, the configuration 0o(x) must have negative 
modes, in other words, it must be unstable. The only candidate for 0o(x) in our case is 
the sphaleron. We arrive at the important conclusion that the real-time instantons describe 
formation of the sphaleron at t — * +oo. It is worth noting that this observation refers, in 
fact, to any ^-instanton which is not real at real time, cf. Ref. [2*%j . 

Let us demonstrate that the real-time instantons determine the suppression exponent for 
the inclusive probability of tunneling from states with given number of incoming particles 
and arbitrary energy 

V(N) = Y,\(f\SP N ^)\ 2 ■ (12) 

U 

From Eq. (jSJ) we obtain, 

V{N) = J d E V(E,N) oc J dE 



e -F(E,N)/g 2 



Let us evaluate the last integral in the saddle point approximation. The saddle point is 
determined by condition 

dE 

On the other hand, formulae (J5j). (fTTjl yield 

9F 

— = -2T . 

dE 

Thus, the integral is saturated by the solution with T = 0, i.e. real-time instanton: 

V(N) oc e~ Frtm/92 . (13) 

The energy E^ rt \N) evaluated on the real-time instanton solution according to Eq. (J8a|) . is 
the best energy for tunneling at given N. We discuss this point further in Sec. 14.41 

Formulae (|12)1. (fT3*j) show that F rt (N) sets the lower bound on the suppression exponent 
F(E) for the one-particle collision-induced tunneling at all energies, extending up to infinity, 
and the smallest value of F(E) coincides with the limit of F rt (N) as N — ► 0. 



10 



2.2 Reformulation of the problem 

In this subsection we adapt the T/#-problem (J7|) for the specifics of the model (JTJ). The key 
observation here is that the characteristic frequency scale of solutions under consideration is 
of the order of the boundary scale /i, which is much greater than m. Thus, the mass term in 
the bulk equation (|7a|) can be neglected, and the general solution in the bulk has the form 

0(t,x) = &(t + x)+0/(t-x), (14) 

where fa and <pf are the incoming and outgoing waves, respectively. They are related by the 
boundary condition (|7b|): 

-0/0) =^sin(&(z) + 0/(z))) ( 15 ) 

where we have promoted t to complex variable z. It is natural to consider 0j and <pf as 
analytic functions of z, and reformulate the rest of the problem (JJJ), i.e. conditions (J7cj) 
and (J7e|) . in the complex z-plane. 

Let us consider the asymptotic future, t — ► +oo (region D in Fig|2j). In this limit the 
field 4>(t,x) is represented by the outgoing wave <pf(t — x) whose argument z = t — x runs 
all the way along the real axis as x changes from to +oo. So, the condition (J7cjl can be 
written in the following way, 

Im (p f (z) = 0, when z<=R. (16) 

On the other hand, it is the incoming wave fa(t + x) which survives in the asymptotic past 
t — > — oo + iT (region A of Fig. EJ), where the argument of the function <pi{z) runs along 
the line Im z = T when x e [0, +oo). So, the condition (|7ej) is reformulated in terms of the 
function 0j. Namely, one performs the Fourier expansion of 0, along the line Im z = T, 

<f>i(z) = J dk Uk)^-^ = J dk {^(ky^-^ +<f )l (-k)e-^ z -^} . (17) 

k>0 

Taking into account that z = t' + iT + x for the initial wave, one compares Eq. (|17|) to 
Eq. (|7Hjl and finds that the positive and negative frequency components and g*_ k of 
solution are proportional to (fii(-k) and <pi(k), k > 0, respectively. Equation (f7e|) takes the 
form 

4>i{—k) = e~ 6 [(j)i(k)]* , A;>0. 
Given this condition, the function fa can be represented as 

fa(z) = X (z - iT) + e~ e ix(z* + iT)}*, (18) 
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where the function 




oo 



dk (f>i(k)e 



ikz 



(19) 



is regular in the upper half plane of its complex argument. Equation (fTSj) provides an alter- 
native formulation of the ^-boundary condition (fTej) . Note that if the number of incoming 
particles is finite, the ingoing wave packet is localized in space. This implies 



[The latter condition is violated in the case of periodic instantons, see Sec. 13.11 ] To sum- 
marize, the T/#-problem is now given by equations ()15j) . (|TH|) and (fT%|) formulated in the 
complex 2-plane. 

To make sure that a solution of the above problem is relevant for tunneling, one should 
check that the value of the field at the boundary x = has correct asymptotics in the 
beginning and the end of the process. Namely, in the case of direct tunneling with the 
soliton in the final state, <p(t, 0) must change from to 2ir as t changes from — oo to +oo. In 
Sec. |Uwe encounter configurations describing creation of the sphaleron at t — > +oo; in that 
case 4>(t, 0) changes from to n. Thus, one obtains the following conditions: 



Finally, let us analyze the issue of existence of a time contour connecting the asymptotic 
regions A and D in Fig. EJ To this end, note that any singularity z\ of the function <pi{z) 
produces a whole half-line of singularities t s = zf — x in the complex time plane, which 
starts from the point z\ and extends to the left parallel to the real time axis, see Fig. |21 
One requires that all these half-lines of singularities of the function <fti(t + x) in the strip 
Im t e [0, T] in complex time plane, should be located to the left of the relevant contour. 
Analogously, one observes that all the singularities of the function <ftf(t — x) in this strip 
should be located to the right of the relevant contour. It is clear that one will always be 
able to find some contour connecting the regions A and D which leaves the singularities of 
(f>i(t + x) and (pf(t — x) to the left and right of it respectively, provided that the singularities 
of these two functions do not coincide. In terms of the variable z this amounts to requiring 
that the singularities, inside the strip Im z e [0,T], of the initial and final wave packets 
<pi(z) and 4>f(z) are situated at different points. Note that this condition is non-trivial, as 
the functions 4>i an d 0/ are related by the differential equation (fTKj) . 



z 



±oo + iT . 



(20) 



4>i(z) + <f>f(z) — > as z — > -co + iT, 
4>i(z) + <j)f(z) — > 2n or it as z — > +oo. 



(21a) 
(21b) 
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Formula (|10j) for the action can also be rewritten in terms of the incoming and outgoing 
waves: 

S = j c dz + <P f M - <P' f ) - - 008(0, + 0/))} . (22) 

The form of the contour C here is somewhat similar to the time contour ABCD in Fig. |21 
it interpolates between the asymptotic regions z — > — oo + iT and 2 — - > +00, leaving the 
singularities of the functions fa and <ftf in the strip Imz e [0,T] to the left and right, 
respectively. 



3 Direct soliton production at low energies 
3.1 Periodic instantons 

To warm up, let us consider periodic instantons which are solutions of the problem with 9 = 0. 
One begins by noticing that Eq. (|15|) possesses a solution, which would be a conventional 
instanton in this model, 



i In 



-i In 



(23) 



The corresponding field configuration 0, Eq. (j!4j) . is real along Euclidean time axis t = —it, it 
changes from to 2tt as r runs from —00 to +00. So, the properties of the configuration (J2*3|) 
are those one expects from an instanton describing vacuum-to-vacuum tunneling. However, 
Euclidean action of the instanton ()23)1 diverges, implying that tunneling at zero energy is 
absent in our case. Obviously, a configuration obtained from ()23|) by the overall change of 
sign is also a solution, which we call ant i- instant on. 

Exact periodic instanton in our model can be constructed with the help of an Ansatz 
inspired by the dilute instanton gas approach [221 123 EI] Let us take a periodic chain of 
alternating instantons and anti- instantons 5 , 
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(24a) 



(24b) 



5 Constants 37r/2 and 7r/2 in Eqs. (|24|l are added for convenience. 
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It is straightforward to check that the functions (f>i, (f)f given by Eqs. ()24|) are real along the 
lines Imz = and Imz = T, so the conditions (|18|) are satisfied with 6 = 0. Note that 
(anti-)instantons in the chain ()24|) are modified by the introduction of a real parameter xq 
which is to be determined from the boundary equation of motion (|15|) . By making use of 
the infinite product representation for the hyperbolic sine (see e.g. Ref. [22], Eq. 1.431.2), 
summation in (|2*4*|) can be performed explicitly: 



3tt 



+ i In 




7T 



— i In 




(25) 



Inserting expressions (|23|) into Eq. ([To]) , we find that the Ansatz satisfies this equation, 
provided that xq satisfies the following relation, 



th 



7TX \ 



1 . 



7T \ T J 

Choosing the contour of integration in expression as prescribed in Sec. l2.2| one evaluates 
the imaginary part of the action: 

2 Im S = 4tt In — + 4tt . 

7T 

Now, it is straightforward to obtain the energy and suppression exponent in the case of 
periodic instanton from Eqs. (JUJ), (fTTjl . 



T 



Att In 



2/t 
~E 



(26) 
(27) 



In accordance with the above discussion, the suppression exponent starts from infinity at 
zero energy. As one could expect, it vanishes at the sphaleron energy E$ = 2/i. 

Let us work out in more detail the structure of the initial and final states of the processes 
described by the periodic instantons. In the initial region (t = iT + t', t' — > — oo) one has 

'n(t' + X + Xq) 



<j)(t', x) = 7r + <f>i(iT + t' + x) = 2 arctg 



cxp 



T 



(28) 



This wave packet has the shape of a kink solution of the sine-Gordon theory with potential 
V oc cos(20). At the same time, the outgoing wave (t — > +oo along the real axis), 

7l(t — X — Xq) 



(j)(t, x) = (j)f(t — x) = 7r + 2 arctg 



exp 



T 
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has the form of antikink. Scattering of such kink-shaped wave packets from the boundary in 
the framework of massless version of the model has been considered in Ref. j23] , and an 
exact expression for the amplitude of kink-to-antikink scattering has been obtained there. 
Our expression (J27|) matches the results of (23] in the weak coupling regime. 

Finally, let us note that in the strictly massless case, the number of particles both in 
the initial and final states of the process described by the periodic instanton is infinite, 
as the field has non-zero asymptotics at spatial infinity, (p(t,x — > +00) = tt. As shown 
in Appendix El introduction of the bulk mass regularizes this divergence and enforces the 
asymptotics (p(t, x — > +00) = 0. The corresponding corrections to the suppression exponent 
(J27j) are small provided the energy of the initial state is much higher than the mass m. 

3.2 Solutions at arbitrary 9 

We now proceed to solutions with 8^0. Let us make the following observation: if (pi is 
real on the real axis, the function <pf determined from Eq. ()15j) with real initial condition, 
is automatically real on the real axis. Thus, all we need is an Ansatz for 0j, which satisfies 
the condition (|18|) and is real on the real axis. One constructs the required Ansatz as an 
appropriate generalization of Eq. ()24a|) . 

+00 

71= — OO 

where xq, Tq are real parameters, with < Tq < T. It is straightforward to check that the 
function 0j determined by ()29|) can be represented in the form (j!8l) . 

Let us show that ()29)) is the most general Ansatz for the function (pi, once the requirement 
of its reality on the real axis is imposed. The function <pf must be regular at the singularity 
point of <pi situated in the strip < Im z < T . This guarantees that the singularity of (pi in 
this strip is logarithmic; we postpone the proof of this statement till Sec. EJ Then, reality 
on the real axis together with the ^-condition (|T5j) fix all other singularities of (pi. Namely, 
due to Eq. (fT%|) any singularity of the function (pi of the form 

C\n(z - z s ) , 

situated below the line Im z = T is accompanied by a singularity 

e- e C* \n(z-2Ti-z*) 



z + x -iT - i2Tn) - i In |(z + x + iT - i2Tn) j , (29) 
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placed symmetrically with respect to the line Imz = T. At the same time, reality on the 
real axis implies that all the singularities of arise in "complex conjugate" pairs situated 
at conjugate points. These two conditions are sufficient to reconstruct the whole chain 
from the logarithmic singularity in the strip < Im z < T. 

The parameter xq in Eq. (J29|) can be removed by a shift of the variable z, and without 
loss of generality we set xq = 0. In order to determine the remaining parameter To of the 
Ansatz, let us analyze equation (fTK|) in the vicinity of the point z = iTq. One represents the 
function 0, in the form 



i In 



[z - lT ) 



+ Ri(z) 



(30) 



where Ri(z) is regular at z — iTq. Power series expansion in Eq. (fTK|) near the point z = iT$ 
has the form 



iTr 



+ R'STo)-<P' f (iT ) + 0(z-tT ) 



-le 



i{Ri{iT )+<t>f(iT )) 



+ M' i (iTa)+i4>' f (iT ) 



+ O(z-iT 



Two leading terms of this equation yield 

e i(Ri(iT )+<l> f (iT )) _ 

R'iiiT ) = 



(31) 
(32) 



These formulae deserve a comment. Considering Eq. ()15|) with given </>j as an ordinary 
differential equation for the function <frf, one might expect all Taylor coefficients of <pf to be 
determined in terms of the function Ri(z) and free integration constant 0/(iTo). However, 
the fact that z = iTo is a singular point of Eq. (fT5j) makes the situation quite different. The 
requirement of regularity of 0/ at this point fixes the value of <pf(iT Q ) according to Eq. (pITj) . 
while the role of the integration constant is played by (p'j-(iT ); besides, the constraint (}3*2"j) 
on the function Ri appears. The latter constraint enables one to determine the parameter 
To, 

T = Ta(6) , (33) 



where the function a(6) is implicitly defined by 

00 „-6»n 

2 « 2 E-r — 2 = 1 - 



n=l 



(34) 
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This relation can be cast into a convenient integral form 6 

[°° sh(qy) 

X e y+ e -i y = ' ( } 

The function 0j is now completely fixed, and 0/ can be obtained by numerical integration 
of Eq. We return to the evaluation of 0/ at the end of this subsection. 

Let us evaluate the imaginary part of the action on the tunneling solution. Surprisingly, 
the detailed knowledge of 0/ is not needed for this purpose. Reality of the solution on 
the real axis implies that the complex conjugate action S* is given by an integral of the 
same function as in Eq. (|22jh but with different contour of integration C*, which is complex 
conjugate to C. Thus, 



2 ImS = -i(S -S*) = -i[ / Cdz - Cdz = -i <t> Cdz . (36) 

\Jc Jc* J Jco 

In the last equality we deformed the sum of the contours C and C* into the contour C Q 
enclosing the singularities 7 z = ±iT of the function 0j. The calculation of the integral (|36| 
is now straightforward: 

2 Im5 = 4vrIm0 / (iT o ) + 4vr . 

Using Eq. (|3*T|) one substitutes Ri(iT ) for cj)f(iT ), takes the former from Eq. (j2Uj) and 
performs summation. The result is 

poo s ^2 fay\ , 

2ImS = 47rln(/zTa)+47r- 16tt / — ■ (37) 

Jo ev+ e - 1 y 

The energy and number of incoming particles are determined from expression (pTFj) in the 
standard way, see Eqs. 0- We find that the energy is given by the same formula (}2l)j) as in 
the periodic instanton case, while the number of incoming particles is 



N = 4tt 



oosh 2 ^ 



dy 



o sh 2 (^) y 



6 One observes that the sum I(a,9) in Eq. I|34() satisfies the differential equation 

d 2 I 2r . 1 



1 



By solving this equation with appropriate boundary conditions, one obtains i|35|) . 

7 In the course of the deformation, the contour does not cross the singularities of (j>f, according to the 
discussion in Sec. 12.21 
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Let us consider in detail the limit 9 — » +00 which is of primary interest. One obtains 

a = 1 - e~ e , (38) 
iV = 47r^e" 9 , -f +00 . (39) 

As expected, the limit of large 9 corresponds to tunneling induced by a few incoming particles. 
The value of the suppression exponent is given by a simple formula in this limit, 

27T/i 



^U=o = 47rln 



E 



(40) 



This is the result we claimed in Introduction for the region I in Fig. ^ 

At first sight, Eq. (piOj) suggests that the suppression vanishes when the energy reaches 
the value 27iyx. In fact, this is not the case: formula (}4T)j) is inapplicable at energies above 
some critical energy E c < 2k /i. The point is that no solution <p/(z) of Eq. ()15j) with required 
properties exists at energies E > E c . Let us clarify this point. 

While the analysis can be carried out for arbitrary 9, it is particularly transparent in 
the case 9 = +00. According to Eq. (jSHJ), the logarithmic singularity of <pi approaches the 
"incoming wave line" Imz = T in this limit 8 , and the Ansatz (|29|) simplifies: 

z — iT 

(t>i — i In — . 

y z + iT 

It is convenient to consider Eq. (fT3j) on the real axis. Introducing u = <fti + <ftf and ( = z/T, 
one writes Eq. ()15|) in the following form, 

— = -A sin w - ^-q— j" > (41) 

where 

A = /iT . 

In new terms, the requirements (J2"T|) . ensuring that the solution is relevant for the soliton 
production, imply the following boundary conditions for u along the real axis: 

u — > 27T , ( — >• -00 , (42a) 
u — > 2,tt or 7r, C +00 , (42b) 



8 This is a generic feature of T/9 instantons: Eqs. (JBJ show that by sending the number of incoming 
particles to zero while keeping their energy fixed, one obtains configurations with infinitely high frequencies 
and thus singular on the initial part of the contour ABCD of Fig|2 We overcome this problem by considering 
the functions 4>i and </>/ in the whole complex plane. 
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where the condition (|42aj) follows from Eq. (j21aj) when one takes into account that the 
asymptotic region z — > — oo + iT and the real axis lie on different sides of the logarithmic 
cut of the function fa. The condition ()42a|) fixes the solution of Eq. (j4~Tj) uniquely. Then, the 
question is whether this solution satisfies the boundary condition (j42b|) . We have analyzed 
this issue both analytically and numerically. Appendix [B] contains an analytic proof of the 
existence of the critical value A c of the parameter A. When A is greater than A c , the solution 
of Eqs. (jUJ), ()42a|) has the correct asymptotics (|42hj) . so the configuration indeed describes 
the production of the soliton. On the other hand, when A is below the critical value A c , 
the asymptotics of the function u at £ — > +00 change to 0. These results are confirmed 
by numerical analysis, see Fig. El As u represents the value of the field at the boundary 




x < x c 



-4 -2 2 4 6 8 10 



Figure 3: Solutions of Eq. (JHJ), found numerically for A = 2.63 (upper curve), A = A c = 2.62 
(middle curve), A = 2.61 (lower curve). 

x = (cf. Eq. (JHJ)), one concludes that the final state of the process at A < A c is the trivial 
vacuum, and no soliton production takes place. At the critical value, A = A c , the function 
u tends to tt when ( — > +00, and the solution describes the production of the sphaleron in 
the final state of the semiclassical evolution. Numerically we find 

A c = 2.62 , (43) 

and the value of the corresponding energy is E c = 1.2Eg- Tunneling at the critical energy is 
still exponentially suppressed: 

F(N = 0, E = E e ) =4vrlnA c . 
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The picture we encounter here has been observed recently in quantum mechanics of two 
degrees of freedom [2E] and in gauge theory [T3j. Results obtained in both cases indicate that 
transitions at energies higher than critical proceed in two stages, creation of the sphaleron 
and its subsequent quantum decay into the relevant final state. The probability of the latter 
process is of order one, while the former is exponentially suppressed. The corresponding 
semiclassical solutions contain the sphaleron at t — > +00. In order to find such solutions one 
has to abandon the requirement of reality of the incoming wave fa on the real axis and thus 
the Ansatz (J2HJ). 



4 Jumps onto sphaleron 

4.1 Solutions with "zero" number of incoming particles 

At energies higher than E c the Ansatz (|29|) is no longer applicable. Still, one can extract some 
information on the properties of (pi from the analysis of Eqs. (fTo'j) . ()18j) . First, let us show 
that the singularity of fa situated inside the strip Imz 6 [0,T] is necessarily logarithmic. 
We denote the position of the singularity by iT . The function 0/ is regular at the point 
z = iT , so <pf and faj can be replaced by constants in small vicinity of this point. Then, 
integration of Eq. (fTK|) yields 



[z) = -f + 2arctg 



1T0 , (44) 



where fo = (pf(iTo), f\ = (p'j(iTo)/ fi, and C is the integration constant fixed by the condition 
that z = iT is the singular point of the function fa. One observes that the singularity of 
the expression (jUj) is logarithmic. 

Thus the function (pi has the form (jHOj) with the function Ri regular in the strip Imz 6 
[0, T]. Besides, one observes that conditions (}3T| . (}32|) on Ri are still valid, as their derivation 
does not make use of any particular Ansatz for fa. Finally, due to Eq. (JTHJ), the presence of 
the logarithmic singularity of fa at z = iT implies the existence of another singularity of fa 
with the structure 



-ie^ In 



2 v z - i(2T - T„)) . (45) 



This information enables one to cast the T/9 problem in the limit 6 — * +00 into the form 
suitable for numerical analysis. As already discussed, this limit is of primary interest as 
it corresponds to semiclassically vanishing number of incoming particles (this means that 
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the real number of incoming particles is much smaller than 1/g 2 )- We now proceed to the 
formulation of the corresponding equations. 

According to Eq. (|18p. the function fa becomes regular in the half-plane Imz > T when 
e~ e tends to zero. As for its singularity in the strip < Im z < T, it hits the line Im z = T in 
the considered limit 9 — > +00, as is clear already from the general reasoning in Footnote |H1 
To demonstrate this explicitly, we write 



i In 



z - iT ) 



- ie~° In 



z-i{2T-T )) +Ri{z), 



where R\(z) is regular both at iTq and i(2T — To). Equation 

-2R t (iT ) . 



then implies 



(46) 



(47) 



T-T 

This formula should be compared to Eqs. (JHHj) . valid for the solutions below the critical 
energy. In that case Ri {iT ) = -1/T ^ 0. Let us assume that Ri {iT ) 7^ also for the 
solutions at energies above E c . Then, Eq. (|4Tj) implies that T indeed approaches T when 
e~ e — > 0. We conclude that in the limit 9 — > +00 the function 64 has the form 



i In 



+ Ri(z) 



(Ah 



where Ri(z) is regular in the upper half-plane. 

As in Sec. 01 it is convenient to consider equations for the functions fa and 0/ on the 
real axis 9 z = x G R. Again, one should be careful about the asymptotics x — ► ±00. The 
function fa satisfies the condition ([20)1 . Taking into account the logarithmic cut of 0« in the 
strip < Imz < T (cf. Sec. 0J) and conditions (|21|). one obtains 



2vr, 
0. 





7T 



as x 
as x 



—00 
+00 



(49a) 
(49b) 



Note that the asymptotics of 0/ at x — +00 corresponds to the formation of the sphaleron 
at the end of the tunneling process. 

A convenient expression for the action functional is obtained in the following way, 

2ImS' = 2Im J dz^fcj)^ — — cos(fa + 0/))) 

/■oo 

= — 47r Im Ri(iT) + 4?r + 2 Im / dx^ffa - - cos(0i + 0/))) . (50) 



3 The reader should not confuse a; which is the real part of the variable z, with the spatial coordinate x. 
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In the second line we used Eq. (|31|) and took the limit To — > T. Integration in the last term is 
performed along the real axis; the first two terms account for the residue of the integrand in 
the logarithmic singularity of fa. Finally, let us present convenient formulae for the energy of 
a solution. One finds two different expressions for the energies of the initial and final states; 
we denote these energies Ei and Ef respectively 10 . For the final energy it is straightforward 
to obtain 

/oo 
dz((j)' f ) 2 , (51) 
-oo 

where the first term is due to the presence of the sphaleron in the final field configuration. 
The analogous expression for the initial energy has the form 

+oo+?.T 

E t = j dz^f. (52) 

-oo+iT 

In the limit e~ e <C 1 the integral in (p)2*j) is saturated by the contribution of the singularity 
at z = i(2T — To). One obtains 

Ei = 2vr-^- = -AttR^T) , (53) 
i — i 

where in the second equality we used Eq. (JJ7J). 

We have performed numerical solution of the following set of equations formulated on the 
real axis: differential equation (JT5j) with the boundary conditions (|49|). condition of reality of 
<pf (Eq. 1)16]) ) and analyticity condition ()48jl . Details of our numerical method are presented 
in Appendix O Let us mention here an interesting property of the solutions. One notes that 
the above set of equations is invariant under the transformation 

fa( z )^2n-(fa(-z*)Y , (54a) 
^(z)^ *-(*/(-**))*. (54b) 

We find that the tunneling solutions at E > E c are symmetric with respect to this transfor- 
mation. This property distinguishes them from the solutions at energies lower than critical. 

The imaginary part of the action, the energy and suppression exponent are calculated 
according to formulae (|5D|). (joTj) . (|53*|) : we use the equality of the initial and final energies as a 
cross-check of the precision of numerical calculations. Results for the suppression exponent 
F(N = 0) are presented in Fig. ^ They cover the interval E c < E < 2.3E$ of the region II. 
10 Evidcntly, Ei equals Ef for any solution of the equations of motion. 
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Let us stress that our numerical solutions describe tunneling with exactly zero semiclassical 
number of incoming particles. The fact that one is able to find such solutions is a peculiarity 
of the model. Unlike in the case of more complicated systems [H] we do not need to perform 
calculations at finite N. 

At energies higher than 2.3E$ the numerical analysis with the method outlined in this 
subsection becomes problematic. The point is that with the growth of energy, the parameter 
T decreases, and the logarithmic singularity of the solution approaches the real axis, accord- 
ing to Eq. (|48j). Hence, the finite-difference approximation used in numerical calculations 
breaks down in the region near z = at small T. Luckily, as we show in the next subsection, 
at \xT <C 1 the singularity can be isolated and treated analytically, thus allowing for complete 
determination of the solution in the limit T — > 0. The energy and suppression exponent of 
the limiting solution yield the point (Ei, F(Ei)) of the graph in Fig. ^ The interval of ener- 
gies between 2.3E$ and E\ is covered by expanding the solution at these energies in powers 
of the small parameter pT around the limiting solution, see Appendix iDl 

4.2 Limit T -> 

In this subsection we determine the tunneling solution in the limit T —>■ 0. We will have to 
resolve the structure of the singularity of the solution in the region near the point z = 0, 
so we keep T small but non-zero, fiT <C 1, and take the limit only in the final formulae 
for the energy and suppression exponent. For the same purpose we introduce small but 
non- vanishing e~ e < 1. It is convenient to choose the contour C, entering the formula (|22|) 
for the action functional, as shown in Fig. EJ From Eqs. (|2()j) . (J2T|) one obtains that the 
asymptotics of the solution on the left part of the contour C is given by 

4>i — > 0, 0/ — > , at x — > — oo , (55) 

while Eq. (|49b|) remains intact. The difference between Eq. (|55j) and Eq. ()49aj) is explained 
by the fact that the contour C and the real axis lie on different sides of the logarithmic cut 
of the function 0j. Correspondingly, the symmetry of solutions considered along the contour 
C is slightly different from (|54|): instead of (|54aJ) one has 

<t>i{z) i ► -<j>*{-z*) . (56) 

We search for solutions invariant under the transformation (j5fij) . (J54bj) . 

Let us introduce e = max{T, e~ e / //}. Our strategy is to separate the complex plane into 
"hard" (\z\ and "soft" (\z\ ^> e) regions, solve equations separately in these regions 
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Figure 4: "Hard" and "soft" regions of the complex plane; respective pieces of solution are 
glued together in the shaded area. Crosses represent the singularities of the function fa in 
the upper half-plane. 



and glue solutions in their intersection e <C \z\ <C see Fig. 0] Let us start with the 
"hard" region. Inserting the expression (|30j) for fa into Eq. (|T5|) and neglecting terms of 
order 0((/xz) 2 ), one obtains 11 

R[-fa f = 1— (l + e WW) . (57) 

Z — 11 o 

An obvious solution to this equation is 

Ri = (f> f = \ . (58) 

The corresponding function 0« satisfies the ^-condition (j!8j) to the zeroth order in e~ 8 . 
However, we saw in Sec. 14.11 that in this approximation the solution is singular on the 
contour C (cf. Eq. (|4~7j) ). So we need to find the regularized solution which satisfies the 
^-conditions to the first order in e~ e . To this end, we linearize Eq. (|57j) : 

r -cp = — . (59) 

z - iT 

11 It is amusing that Eq. (|57|l is symmetric with respect to conformal transformations 
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Here notations r = R4 — n/2, p> = 0/ — vr/2 are introduced. Equation (|59j) can be solved 
with respect to any of the two functions, ip or r. One gets 

2r'(z 1 ) 



<p — — r + (z — iT 
2 

r = (p — 



JiT 



Z\ - iT 
(p(z 1 )dz 1 , 



dz\ , 



(60) 
(61) 



where the integration constant in the expression (JfiT|) is chosen in such a way that the function 
r is regular at the point z = iT . We do not specify the value of the integration constant in 
Eq. ()60|) as this constant is irrelevant for what follows. 

Let us investigate the singularities of the functions <p, r. In the same way as in the 
previous subsection, one observes that the only singularity of r in the upper half-plane is 
situated at z = i(2T — T ) and has the form (|3Hj) (cf. Eq. (JIHJ)). Then, it follows from 
Eq. ()6()j) that the only singularity of 'p(z) in the upper half-plane is situated at the same 
point, z = i{2T — T ) (an apparent singularity at z — iT is eliminated, as before, by 
imposing condition ()H2j) on the function r). Reality of ip on the real axis implies now that 
the only other singularity of p> is situated at the point z = —i{2T — T ), its structure being 
"complex conjugate" to that of the singularity at z — i(2T — T ). We determine the form of 
these singularities from Eq. ()6()j) . single them out explicitly and expand the remaining entire 
function into power series. In this way we get the following representation: 



<P = - 



_ y fi(z - iT) In [|(* - i(2T - T ))] + p{z + iT) In [^(z + i(2T - T )) 
7rTe~ 



T-Tr 



+ A fxz + ^A n {^z) 



l+2n 



(62) 



n=l 



The first term in the second line as well as the presence of only odd powers in the power 
series are fixed by the symmetry (|54bjl . The coefficients in Eq. (}6*2*|) are to be determined 
by matching this expression to the solution in the "soft" region. We will see below that the 
function 0/ in the latter region is of order one, and thus all the coefficients, ^t-t ) > ^ n 
are of order one. So, the contribution of the sum in Eq. (|52Jl can be neglected in the "hard" 
domain; we omit it in what follows. Inserting expression ()62j) into Eq. ()61jl one gets: 
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T-Tn 



i{2T-T )) 



, T + T 
+ 1— — — e 



Z - zT n 



■ 111 



1 + 



T- 

iT 



T 



In 



z + i{2T-T Q )) 



2iT 



T-Tn 



(z + tT ) + 



7rTe- e 

T-Tn 



(63) 



iA fiT . 
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In the region e <C \z\ <C the expressions (|62p. ()63j) simplify. In terms of the original 
functions </>j, </>f, one obtains 



z In 



LI ~\ 7T C u 7T £ u I LI Z 2 

-(z-iO) H — — —z, 6*= — — — zln - — ] + A Q uz , (64) 

2 V '\ 2 T-T ' r/ 2 T-T V 4 / 

where terms of order O(fiT), 0(e~ e ), 0((/iz) 2 ) have been neglected. 

We now turn to the "soft" region. Here the contour C coincides with the real axis. 
Introducing the real and imaginary parts of the incoming wave, (f>i(x) = a(x) + ib(x), one 
rewrites Eq. (|THj) as a set of two real equations, 

b' = fishbcosu , (65) 
u' = 2a' — /j, ch 6 sin u , (66) 

where u(x) = a(x) + 0/(x). In the "soft" region one can neglect the regularizing parameters 
/iT, e~ e . [Solutions at non-zero T are considered in Appendix [DJ] Then, the function (pi 
is regular in the upper half-plane of the "soft" region, and the derivatives of its real and 
imaginary parts are related by an analog of the Cauchy formula, 

a'(x) = — V.P. f°° d Xl Xlb '^ , (67) 

71X J —oo ^ 

where the integral is understood in the sense of principal value. Equations (JHHJ), (|B7jl 
are supplemented by boundary conditions at x = which originate from matching with the 
"hard" part of solution. From the expressions (J64|) one obtains 

b — > In (t^ x j j a_ ^^"; u ^ rr , at x ^ +0 . (68) 

The conditions at x — > —0 can be reconstructed from the requirement of symmetry with 
respect to ()56|). (|54bj) . For completeness we write down the boundary conditions (|49b[) . f|55jl 
explicitly: 

a ^ , 6 — > , it — > tt , at x — > +oo , (69) 
while the symmetry transformations ()54b|) . (j3B]) are: 

&(x) i— > £>(— x) , a(x) I— >■ — a(— x) , w(x) i — > 7r — tt(— x) . (70) 

Due to the symmetry (|7()jl . it is sufficient to consider half of the real axis, and we assume 
x > in what follows. Given the function u, equation (|65|) can be integrated explicitly: 

b(x) = lnth f — ^ / cos«(xi)dxi J , (71) 
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where the condition (|68|) is imposed on b(x) at the origin 12 . Equations 1)71)) . (|67|). ()66|) form 
a complete system. 

One obtains a fairly good approximation by linearizing Eqs. f!71|) . with respect to 
5m = u — tt. Equation (J71|) takes the form 

b = lnth ( —x 



Substitution of this expression into Eq. (jHIJ) yields 



■*(>( 



1 



TT 



V 7T 



where 
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x + 1 



In 



1 'f 



Now, <5m can be found from Eq. ()66|). The corresponding analytic expression is not illumi- 
nating, and we do not present it here. The linear approximation is improved by numerical 
iterations. Each cycle consists of solving Eqs. (J7T|) . (pTTjl . (jfifij) thus determining functions 
6, a, 5m one after another, starting from the approximation for 8u obtained in the previous 
cycle. After 30 iterations one obtains numerical solution with precision of order 10~ 6 . The 
resulting functions a, b, 5u together with those obtained in the linear approximation are 
plotted in Fig. 

Let us consider the behavior of the solution at x — > +0. Equation (jfifij) allows for the 
following behavior: 
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'- - A/j,x + 0(x 2 ) 



- - 2A{ix In (J^-J + B/ix + 0[ 



x 



(72) 



We observe that both the solution in the linear approximation and the full numerical solution 
exhibit this behaviour. Numerically, A = 0.845, B = —2.90. Matching these expressions to 
the asymptotics (|64|) of the "hard" part of the solution, one obtains 



T-T n 



Afj, , A = B . 



(73) 



Thus, the limiting solution is completely determined. 

A comment is in order. Strictly speaking, in the limit T, e~ 9 — ► the "soft" region covers 

the whole complex plane and one is left with the "soft" part of the solution only. However, 

12 Note that due to the asymptotics 1)69(1 of u(x) one has b oc e~^ x at large values of x. Such behavior 
indicates that at t — * +oo the imaginary part of the solution describes the evolution along the negative mode 
of the sphalcron. This is a generic feature of tunneling solutions describing jumps onto the sphaleron |28|. 
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Figure 5: Limiting solution in the "soft" region. Dashed lines represent linear approximation 
in 5u. 

the solution is then singular at z — 0, and it is impossible to calculate its energy and action. 
So the term "limiting solution" implies not only the form of the solution in the "soft" region 
but also the appropriate resolution of the singularity, Eqs. (JS2J), Ijfi3jl. 

The next step is to evaluate the imaginary part of the action functional on the constructed 
solution. It is natural to separate it into the sum of contributions of the "hard" and "soft" 
regions. Neglecting terms of order fj,T, e~ e we have 

in 1 



2lvaS hard = 21m 



A* 



dz = 2n . 



C,\z\<x, 



2{z — iT ) z — iT 

Here is some number in the range e Ci» C For the "soft" contribution we write 
2lmS so f t = 2 / {(j)fb' — /ishfrsinn) dx 

J |x|>x* 

POO POO 

= 2nb(x*)+4 (cj) f — tgu)tf dx = 4 
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cos 2 u 



(74) 



0j J b dx , 



•x* JO 

where in passing from the first to the second line we used Eq. (J65j) and symmetry with respect 
to (|7L)j). while in the last expression terms vanishing in the limit x* — > are neglected. A 
good estimate for 2 Im S so f t is obtained in the linear approximation: 

In 2 lnvr' 



2 lmSf? ft = 4 J a'bdx = 4vr ^— 6C'(— 1) 



3.83 



6 2 

where integration is performed using Eqs. 8.371.2, 3.951.14 from Ref. [321- The calculation 
of the integral (|71jl for the full numerical solution yields 2lmS so f t = 3.98. Note that this 
value is comparable to the contribution 2tt coming from the "hard" core. 
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An important quantity is the energy of the limiting solution. Again, we calculate sepa- 
rately the energies of the initial and final states. While formula (|51|) for the final energy is 
directly applicable in the limit T — > 0, the analog of Eq. (p)3|) for the initial energy is 

Ei = 2tt— — = 2nAfx , 
l — i 

where A is defined in Eqs. (I72J) . For our numerical solution Ef = Ei with good accuracy, this 
may be viewed as a confirmation of the consistency of the numerical method; the obtained 
value is 

E l = 2.65E S . (75) 

At first sight, the result (J75j) is puzzling. Indeed, one could expect to cover the whole range 
of energies up to infinity by sending the parameter T to zero. However from Eq. (|75j) we 
learn that this is not the case: T/#-problem does not have any solutions with energies higher 
than Ei. On the other hand, the suppression exponent is still non- vanishing at the limiting 

energy, 

F(N = 0, E = Ei) = 10.27 . (76) 

A question is then, how tunneling occurs at energies above Ei, and what is the corresponding 
tunneling probability. 

The answer is suggested by an observation that in all the formulae describing the limiting 
solution, even in those related to the "hard" core, one can set T = (it is sufficient to keep 
e~ e non- zero to resolve the singularity of the solution). As discussed at the end of Sec. 12.11 
the solution obtained in this way belongs to the class of real-time instantons. Thus, it gives 
the maximum probability for tunneling from states with semiclassically vanishing number 
of particles. Consequently, the suppression exponent F(N = 0, E > E{) is not smaller than 
F(N = 0, E = Ei). It cannot be larger either, because one can always imagine a tunneling 
process which has the same exponential suppression as the process at the limiting energy. 
Namely, incoming particle(s) can release the energy excess {E — E{) by the perturbative 
emission of a few other particles at the beginning of the process, so that subsequent tunneling 
occurs effectively at the energy E\. From this physical reasoning one concludes that the 
suppression exponent stays constant at energies higher than Ei (cf. Ref. |15j). Remarkably, 
the semiclassical treatment is still possible at these energies. We return to this issue in 
Sec. 14.41 Now, let us turn to the real-time instantons which, for their novelty, deserve a 
somewhat detailed study. 
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4.3 Real-time instantons 



Let us recall that the real-time instantons are solutions of the T/^-problem with T = and 
finite 9. In the model under consideration these solutions can be obtained by an iterative 
method which is a slight modification of that used in the previous subsection. As before, 
one introduces the functions a, b, u satisfying Eqs. Note that, unlike in the case 

of Sec. 14.21 all the functions are supposed to be regular on the entire real axis, including the 
point x = 0. Requirements of regularity together with symmetry under the transformations 




Figure 6: Example of real-time instanton configuration with 6 = 2.04, N = 2.3. 

(jTHj) lead to the following conditions at x — 0, 

6'(0) = , a(0) = , u(0) = vr/2 . (77) 

Formula (Jfi7|) is not directly applicable in the case of finite 9. An appropriate modification 
comes about when one notices that a direct analog of Eq. (joTj) relates the real a and imaginary 
b parts of the function x introduced in Eqs. (fTSjl. (JTHJ). Using Eq. (fTHj). one expresses the 
functions a, b in terms of a, b and gets, instead of Eq. (fH7]l . 

a (^J = T ' — V.P. / axi . (78) 

1 — e y 7ra; y_ oc X\ — x 

Finally, Eq. (|65p can be integrated explicitly in the present well, but now the inte- 

gration constant is arbitrary, 



b(x) = lnth ^— ^ J cosu(xi)dxi + x 



(79) 
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Figure 7: Energy and suppression exponent versus the number of incoming particles in the 
case of real-time instantons. 

We obtain numerical solutions by iterating Eqs. (ff9|) . (fT8j) . (}66|) in the indicated order. The 
value of x is determined at each cycle of iterations from the requirement that the function u 
satisfies the boundary conditions (|77jl . (jfi9j) . An example of real-time instanton configuration 
is shown in Fig. |U1 

It is straightforward to obtain the energy and imaginary part of the action from Eqs. (Jo"Tj) 
and ()22j) . As for the number of incoming particles, in the case of real-time instantons it can 
be expressed in a particularly convenient form, 



which can be obtained from Eq. (j8b|) in a straightforward way. The energy and suppression 
exponent as functions of the initial number of particles are shown in Fig. The curves start 
at iV = from the values corresponding to the limiting solution and end up at iV = 3.7 where 
the suppression exponent is zero. We see that the branch of real-time instantons interpolates 
between the limiting solution and the region of classically allowed transitions. 

4.4 Tunneling above the limiting energy 

In this subsection we follow an approach which is somewhat different from that used in the 
rest of the paper. Instead of considering the inclusive probability (JHJ), we study tunneling 
from a single initial state with given number of particles N and energy E. The energy is 
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assumed to be higher than the energy of the real-time instanton corresponding to the same 
initial number of particles, E > E^ rt \N). Let us demonstrate that the suppression exponent 
of tunneling from an appropriately chosen initial state coincides with that of the real-time 
instanton. 

Using the approach of Ref . [231 , one can show that any field configuration which satisfies 
Eqs. (|7a| . (J7EJ), (J7cj) (but not necessarily Eq. (|7ejl ) in the complex time plane determines the 
suppression exponent for the tunneling probability of the form 

?(£,iV; fl ) = V^^. (80) 
j (a\P E PN\a) 

Here summation is performed over the final states only, the initial state PePn\o) is a coherent 
state projected onto the subspace with fixed energy and number of particles. The coherent 
state corresponding to a given solution is defined by 

Ofc|a) = f k \a) , (81) 

where /& are the positive frequency components of the solution x) at t — > — oo, 

1 f dk 



x 



/dk 
—== (f k e- iUkt+ikx + g* k e iuJkt - ikx ) . 



Note that unlike in Eq. (|7d[) we now consider the asymptotics along the real t-axis. The 
suppression exponent for the probability (jSUj) is still given by Eq. (fTT|) where the parameters 
T and 9 are now determined from the following equations, 



E = J dku k f k fZe^ 1+ » , N = J dkf k r k e^ 1+y . (82) 

One can check that, as before, E coincides with the actual energy of the solution, and the 
expressions (JHJ for the energy and number of incoming particles are applicable. Note that 
in general the average energy and number of particles in the state \a) are different from E 
and N. 

We consider the following configurations which differ from the real-time instanton by a 
small perturbation: 

M*) = ^/ + (83) 

where the perturbation is represented by the function ip(z) which is real on the real axis 
and falls off fast enough at x —>■ ±oo. When the parameters a, j3 are small, ^ < a < 1, 
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the configuration (|83j) is an approximate solution of Eq. (|15|). Evaluation of the energy and 
number of incoming particles using Eqs. (|51|). (|8b|) yields 



the solution ()83|) describes tunneling from a state which differs from the initial state of the 
real-time instanton by addition of a few particles carrying a finite amount of energy. Compar- 
ison of the structure of the initial and final waves shows that these particles scatter elastically 
from the boundary x = and do not affect the tunneling process. It is straightforward to 
evaluate the imaginary part of the action of the solution (|H3*|) and determine the parameters 
T and 6 from the formulae (|82jl . In this way one finds that the suppression exponent of the 
solution ()83|) coincides with that of the real-time instanton up to terms 0(a). 

As before, the suppression exponent of collision-induced tunneling is obtained in the 
limit N —>■ 0. One concludes that it is equal to that of the limiting solution at all energies 
above E\. 

5 Concluding remarks 

The results obtained in this paper show that in the model (0) the soliton production due to 
collision of a particle with the boundary is always exponentially suppressed. This may be 
viewed as an evidence in favour of the exponential suppression of collision-induced tunneling 
at all energies in other theories. We have identified the mechanism of suppression at high 
energies of colliding particles. It is a consequence of substantial rearrangement the system 
has to undergo in the course of transition. Moreover, this rearrangement occurs in an optimal 
way at a certain limiting energy Ei where the suppression is the weakest. At higher energies 
the system emits the excess of energy in the form of particles which do not interfere with the 
tunneling process, and subsequent transition occurs with the energy E t . The consequence 
of the latter phenomenon is that the suppression exponent stays constant all the way from 
E = Ei to E = oo. It is natural to conjecture that this mechanism of suppression is generic 
for collision-induced tunneling in quantum field theory. 



E = E {rt) + 




N = N {rt) + 0{a 2 ) . 



These formulae demonstrate that in the limit 



ot,(3— >0, a 2 / (3 fixed 



33 



As a by-product we have found that there is a novel type of tunneling configurations which 
we call real-time instantons. These are solutions to the T/#-problem with the parameter T 
set equal to zero. From the physical point of view, real-time instantons provide maximum 
probability for tunneling from initial states with given number N of incoming particles. In 
the plane of parameters (E,N), the line of real-time instantons interpolates between the 
region of classically allowed transitions and the point (E = Ei,N = 0). The suppression 
exponent calculated for a real-time instanton F rt (N) sets a lower bound on the suppression 
exponent F(E) for the tunneling induced by one or a few incoming particles at all energies. 
Modulo standard assumptions, F rt (N) in the limit N — > coincides with F(E) at energies 
higher than E\. Thus, one can use real-time instantons to determine the value of the limiting 
energy E\ and the suppression exponent of collision-induced tunneling at energies higher than 
Ei in various models. 

Finally, let us recall that the model investigated in this paper describes various condensed 
matter systems. It would be interesting to work out implications of our results for these 
systems. 
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A Periodic instantons at non-zero bulk mass 

Let us consider the effect of non-zero bulk mass m on the periodic instantons found in 
Sec. 13.11 We assume that the mass is small, mT -C 1. An approximate periodic instanton 
4> m (t,x) in the massive case is constructed in the following way. In the region \t\ <C 1/m, 
one has 

m (t,x) = o (t,x)e- mx , (84) 
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where (fio(t,x) is the periodic instanton solution in the massless case given by Eqs. (|T1|). (}2l)j) . 
It is straightforward to check that the configuration (J84j) satisfies the field equations (J7a|) . 
(|7b|) up to terms of order 0(mT). At \t\ ^> T the value of the field at the boundary is 
frozen to and 2tt in the parts A and D of the contour ABCD, respectively (see Fig. |2J). 
Let us concentrate on the initial part of the contour; we set t = t' + %T . In the region 
T <C \t'\ <C 1/m, the configuration ()84j) should be matched to the solution of the free 
massive equation, which has the form ()7djl . Making use of Eq. (|28jl . one obtains the Fourier 
transform of the function (f) m : 

f 00 rib _,p ifa o , i 

<U*\ x) = / • -^re*(* , T « |f | « 1 



«m 2 cn m 



2ch(f) 

where we omitted terms negligible in the limit m — > 0. Comparing this expression to Eq. (|73l 
we get, 



fk = 0(-k) 

gt = e(-k) 



(k - im) ch (^f) 
{—k - im) ch (^) 



Now, it is straightforward to evaluate the energy and the number of incoming particles using 
Eqs. (jSj) . The energy is given by the same formula as in the massless case, Eq. (|26j) . while 
the initial number of particles is 

f°° 7T 

WTWaW) dk -* HmT)+0(1) - 

We see that the small bulk mass does not affect the energy of the periodic instanton. On the 
other hand, it regularizes the number of incoming particles which diverges logarithmically 
in the limit m — > 0. To find the effect of the mass on the value of the action, one performs 
integration by parts in expression (fTU|) and obtains 

S = J dt (^-(pd x (p — — cos 4>)^j 

From the above consideration it follows that the integrand is non-vanishing only at t ~ T 
1/m. So, one substitutes the expression (j84j) into this formula and finds that the value of 
the action differs from that in the massles case by terms of order 0(mT). We conclude that 
the influence of the small mass on the suppression exponent of tunneling described by the 
periodic instanton is small when the energy of the process is much higher than the mass. 
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B Analytic proof of existence of the critical energy 



We are going to prove the following statements concerning Eq. (|41|): 

a) if A > 4, a solution of Eq. (J4*T|) with asymptotics (J4~2|) exists; 

b) if A < 7r/2, there is no solution of Eq. (j4*Tj) with asymptotics (|32~|l . 

We start with the statement (a). At A > 4 the lines defined by equation 

4 

- - ^ + A sin u = (85) 

separate the band tt < u < 2tt of the (C,u) plane into three regions, where the derivative 
of solution of Eq. (J41j) has different signs, see FigJHJ The condition (J42a|) picks up a single 

u 



2n 




c 



Figure 8: Regions of positive and negative u' and qualitative form of the solution of Eq.([4ip 
in the case A > 4. 

solution of Eq. (j41j) . The integral curve of this solution starts from the region I at ( = — oo. 
As £ increases, it crosses the line u' = and comes into the region II where u' > 0. Once 
in the region II, the integral curve cannot leave it. Thus, u' stays positive all the way to 
( = +oo, and the integral curve gets attracted to the asymptotics 27r when ( — > +oo. 

We now proceed to the proof of the statement (b). At A < 7r/2, the form of the lines 
(|83|) is different, see Fig. El The tips of these lines are situated at ( = ±Co, where 

Co = ^4/A - 1 . 

At ( < — Co the integral curve of the solution lies in the region I where u' < 0. Let us show 

36 




Figure 9: Regions of positive and negative u' and qualitative form of the solution of Eq.([4ip 
in the case A < 7r/2. 

that it crosses the line u = tt from above within the interval — Co < C < Co- Indeed, Eq. (|41j) 
implies, 

> ^ ^ 4 
Integration of the latter inequality yields 

w(Co) < "(-Co) + 2ACo - 8 arctgCo < vr . 

Once the integral curve comes below the line u — tt, it cannot reach the region II where 
u' > (see Fig. EJ). Thus it does not approach the asymptotics 2n at ( — > +oo. 

The above analysis necessitates the existence of the critical value tt/2 < A c < 4 at 
which the behavior of the solution changes qualitatively. We have obtained the value of A c 
numerically, see Eq. (|4"3~j). 

C Numerical method at E > E c 

In this Appendix we formulate a discretized version of the problem (fT5|). (fTBj). (|4"5|). P9*)l. 

When searching for an appropriate numerical method, one encounters two major difficulties. 

First, since all the relevant solutions at E > E c are close to the sphaleron configuration at 

t — > +oo, they display an instability similar to that of the sphaleron 13 . To overcome this 

13 Instability of this kind is a generic feature of tunneling solutions at energies higher than critical, see |2H1 
El El- A general method of dealing with this instability is proposed in |28j . 
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difficulty, we use specifics of the system in hand. As already mentioned, equations (|15|) . 

(J48j) and boundary conditions (|49|) are invariant under the transformation (J54|) . So, 
it is natural to look for solutions symmetric with respect to this transformation. Once the 
symmetry requirement is imposed, the set of considered configurations is restricted to those 
close to the sphaleron at t — > +00, and solutions become stable. 

The second numerical challenge is the slow falloff of the solutions at x —>■ ±00. Equation 
(|15j) allows for solutions with asymptotic behavior 1/x at x — > ±00, and numerical analysis 
confirms that this is indeed the case. We overcome this problem by introducing a non- 
uniform lattice with link size growing when x — > ±00. This is achieved by performing 
conformal transformation of the complex z-plane, 

z — iT , . 

w = = , 86 

z + iT ' v ; 

which maps the real axis onto the the unit circle w = e ta . In particular, the asymptotic 
regions x — > —00 and x — > +00 correspond to the vicinities of the points a = and a = 2tt, 
respectively. It is clear that the uniform lattice in angle a produces discretization of the 
required type of the real z-axis. 

We proceed by reformulating the problem on the circle = 1. First, it is convenient to 
introduce the function 

(z iT 
z + iT 

Due to Eq. (|48|h S{(z) is regular in the upper z-plane. Thus, after the conformal transfor- 
mation (J86j) it is regular inside the unit circle \w\ < 1, and the Taylor series 

00 



Si(w) = i^^SkW* 



k=0 

converges everywhere inside this circle. For w = e lot we obtain 

00 

S l (a)=iJ2Ske lka . (87a) 

k=0 

The requirement that the Fourier series of the function Si has the form ()87a|) is equivalent 
to the condition (J48|) . Equation (|15|) in new terms reads 

^[ S --^= 1+ 2i^ Sm(S ' + ^- a) - (8?b) 

where the function 0/ is real due to Eq. (jlfij) . The reformulation of the boundary conditions 
(|49|) is straightforward. However, all but one of the conditions (J4"9"j) are redundant. Indeed, 
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if the conditions at x — > — oo (a — > 0) are satisfied, requirements at x — > +00 (a — > 27r) 
hold automatically due to the symmetry (|54|) . Moreover, the condition on the function Si 
at a — ► is satisfied automatically due to Eq. (|87b|) . The remaining condition implies 

</>/^0 as a -> . (87c) 

One notes that the symmetry (jo^j) relates the points a and 2ir — a of the circle |w| = 1. 
Therefore, Eq. (jH2h|) can be regarded as a differential equation on the orbifold a G [0, 7r], with 
the condition (|87c|) imposed at the orbifold boundary. One also notes that the transformation 
(|54|) brings the Fourier image Sk into its complex conjugate S^, so that the invar iance of the 
function Si under the symmetry (|SH) implies that G JR.. 

To discretize the problem ([87)1 . one introduces the uniform lattice with sites at the points 
a n = A(n + 1/2), where A = tt/N, n = 0, . . . , N — 1. The boundary condition (|87cj) takes 
the form 

(<P f ) n=Q = . (88a) 

The condition at the other boundary of the orbifold, a = n, is obtained in the following way. 
One discretizes Eq. (18 7b J) in a straightforward manner and writes down a finite-difference 
equation relating (N — l)-th and iV-th sites of the circle. Then, the values of the fields at 
the iV-th site are expressed in terms of the values at (N — l)-th site by using the orb ifo Id- 
symmetry condition, and one obtains 

^ [tt - 2(0/^-1 + 2 Re^iv-i] = 

1 + o ■ 2, \ sin ( Re(Si)jv-i + ((fif)N-i - «jv-i) cosh ( lm(S i ) N - 1 ) . (88b) 

2 sin (ajv_i/ 2 ) 



When discretizing Eq. ()87b|) in the internal points of the orbifold, one has to deal with the 
explicit singularity at a = in the right hand side of this equation. After several attempts, 
we have found a discretization scheme which does not produce numerical artifacts at a = 0. 
It has the following form 

^[(Si) n - (5i) n _i] - ^ [(#/)„ - (0/)n-l] = 
1 + o • 27 777: \ Sil1 (('S'On-l + (^/)« ~ a n-l/2) + COS ((S^n-i + (0/) n - a„_l/ 2 ) 

I sin [a n -i/2/ 2)1 * 
/i 2 T 2 A 2 

sin ((5i) n + (0/) n - a n _i/ 2 ) cos ((-Sijn-a + (<f>f) n - a n _i/ 2 ) • (88c) 



4sin 4 (a n _ 1/2 /2) 
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Here the index n runs from 1 to iV — 1. It is straightfoward to check that Eq. ()88c|) for a n ~ 1 
provides a second-order approximation of equation (j87b|) . At the same time, for a n ~ A the 
last term of Eq. ()88c|) becomes dominant, and one obtains 

(Si) n + (<f) f ) n -a n = when a n ~ A . 

Again, the discrepancy between this equation and the continuous one is 0(A 2 ). The approx- 
imation turns out to be worse in the region between a n ~ A and a n ~ 1. Careful numerical 
analysis of the scaling of error with the lattice spacing shows that the approximation order 
lies somewhere in between 1 and 2. 

Finally, the Fourier transform (|87a|) is replaced by its discrete version, with the sum 
taken from k = to k = N — 1. We choose the real variables Sk and (0/) n as the unknowns. 
Since k, n = . . . (N — 1), the total number of unknowns, 2N, coincides with the number of 
real equations in the system (|K%|) (two real boundary conditions ()88a|) . (J88bj) and (N — 1) 
complex equations (j88c|0 . We solve the system of equations (JHHj) by the Newton-Raphson 
method. Then, the values of the final energy and imaginary part of the action are given 
by the formulae obtained by a straighforward discretization of Eqs. (|5U|) . (|51|). In order to 
calculate numerically the initial energy (p)3j) . we note that 

B&z = iT) = -L [1 - iS'^w = 0)] = ± [1 + SW] • 

The comparison between the values of the initial and final energies yields the estimate of 
the discretization error. In our calculations we used the lattice with 2048 sites, and the 
discretization error was smaller than 10~ 3 for all our solutions. 

D Solutions at small T 

The purpose of this Appendix is twofold. First, we are going to clarify the structure of the 
solution of the T/#-problem at small values of the parameter T. Second, we will obtain the 
formulae for the suppression exponent in the region of energies 2.3E$ < E < E\ which is not 
covered by the numerical method of Sec. 14.11 To make analysis more transparent, we work 
below in the limit e~ e = and set \i = 1. 

First we observe that the limiting solution found in Sec. 14.21 represents the first terms 
of a systematic expansion of the solution in powers of the parameter T. We start with the 
"hard" core. One notices that the characteristic length scale in the "hard" region is set by 
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the parameter T. So, it is natural to represent the solution in this region in the following 
form, 

€f = </ 0) (0 + (C) + T 2 <^ 2) (C) + ■ • • , (89) 

where ( is a rescaled variable, 

C = z/T. 

We allow the functions (f)[ j (C)j <t>\ /^(C)> • • • to contain terms of order (lnT) n . To illustrate 
this point let us write the "hard" core solution of Sec. 14.21 in terms of the rescaled variable (. 
Making use of Eqs. d), one obtains, 



i + zln (l) + zln (c-^) 

T{2zAlnT + 2zA^^ln ) - A( - i(A + B) \ (90a) 



(90b) 



(-i \ 2i 

dh) * ~ f ■ ^ - ( T 



f — — + T < -A(C - i) ln(C -%)- A(C + i) ln(C + i) - 2A( In I - ) + B( - nA 



where the coefficients A and B are defined in Eqs. (j72j) . It is clear that these expressions 
represent the first two terms of the expansion ([89)1 . 

In the "soft" region the characteristic length scale is of order 1. So, the expansion in the 
"soft" region has the form 

4® = 4f(z) + T$}\z) + T*4f(z) + .... (91) 

The first term of this expansion has been found in Sec. 14.21 

To construct the solution in the entire complex plane one should match the asymptotics 
of the expansion (|KH|) at £ — > oo with the Taylor series of the expansion ()91)1 at z — > 0. Note 
that because the arguments in the functional expansions (jHHJ) and (|9*T|) differ by rescaling, 
matching intertwines terms of different order in T in the expansions (|89jh ()9ip. Thus, in 
Sec. 14.21 we have used the behavior of the zeroth order approximation in the "soft" region 
to determine the coefficients of the "hard" part of the solution to the first order in T, see 
Eqs. f!73|) . Here we proceed to calculate the first order term in the "soft" part of the solution. 

Since e~ 9 = 0, the function (f)f 1 ^ is regular in the upper half-plane of the "soft" region. 
Thus, its real and imaginary parts on the real axis, a^'(x) and b^(x), are related by the 
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Cauchy formula (|67|). Equations (|65p. ()66|) in the first order approximation are linearized: 

&(*) = -s(x)u {1 \ (92) 

(93) 



d , ' 
dx- ~ C(X) 



dx 



+ c x 



"> 2 — - S (o;)&« 



d.T 



Here, as usual, = + while the coefficients c(x), s(x) are 



clx) 



ch&<°> 



cosu 



(0) 



six) 



sin w 



(o) 



Boundary conditions at x = originate from matching with the "hard" part of the solution. 
One writes expressions ()90|) in terms of the variable z = T( and extracts the corrections pro- 
portional to T, by formally considering the variable z as a quantity of order one. According 

(si) 

to the above discussion, this fixes the small-z behavior of the functions 4>\ f ■ 



{i + 2zAln(| + iO) -i{A + B) + nA} 



-TlA 



at 2 







Finally, we recall that due to Eqs. ()49|) the finctions vS 1 ' vanish at x — ■> +oo. 

The problem outlined above can be solved by iterative method similar to that of Sec. 14.21 
Namely, one starts from «W = 1/x and finds the functions by solving 

Eqs. f|92|) . (|67|). (J93|) one after another. In this way the improved approximation for 
is obtained, and the next cycle of iterations begins. After 30 iterations the solution is found 
numerically with accuracy of order 10~ 4 . 

We have determined the solution in the first-order approximation both in the "soft" and 
"hard" regions. So, we are able to calculate its energy by performing the integration in 
Eq. ()52|) in the first-order approximation. We find 



E(T) =Ei + Air AT (2A\nT - B) +WT . 



(94) 



In this expression E\ stands for the limiting value, the second term is the contribution of 
the first-order corrections in the "hard" core, and the last term represents the correction 
calculated numerically in the "soft" region, 



W 



dx 



9.6 
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Given the expression (|94|) for the energy, it is straightforward to evaluate the suppression 
exponent using Eqs. ©, 



F| JV=0 = F / -r 2 [/ 1 lnT + / 2 ] J (95) 

where Fi is the suppression at the limiting energy; the constants fi, f 2 are related to the 
coefficients of Eq. flU)). Numerically, fx = 17.9, f 2 = 49.5. 

We use the expressions ffHljl . (fHKj) to determine the function in the region 

2.3Es < E < Ei. Note, though, that the expressions (}94j) . (}95|) are valid only in the vicinity 
of the point E = Ei; corrections to these formulae grow when one goes away from the lim- 
iting energy. At E = 2.3E$ we are able to compare the value of the suppression exponent 
given by Eqs. (jHlJ), (|93jl with the one calculated numerically by the method of Sec. 14.11 We 
have found that these two values coincide with the precision of order 0.5%. 
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